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We present the first results in the search for relativistic magnetic monopoles with the IceCube 
detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume 
of 1 km 3 . This analysis searches data taken on the partially completed detector during 2007 when 
roug hly 0.2 km 3 of ice was instrumented. The lack of candidate events leads to an upper limit on 
the flux of relativistic magnetic monopoles of $go%cx. ~ 3 x 10~ 18 cm~ 2 sr -1 s -1 for /3 > 0.8. This 
is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost 
7 below 10 7 . This result is then interpreted for a wide range of mass and kinetic energy values. 

PACS numbers: 14.80.Hv 



I. INTRODUCTION 



Magnetic monopoles are an important element in a 
complete picture of our universe. Their existence would 
explain the quantization of electric (and magnetic) charge 
via the Dirac quantization equation g — Ne/2a [lj. 
They appear as topological defects from symmetry break- 
ing m Grand Unified Theories (GUTs) 0] with masses 
~ 10 4 — 10 17 GeV [3|, depending on the breaking scheme. 
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Additionally, they would bring a complete symmetry to 
Maxwell's equations. 

Magnetic monopoles produced in the early universe via 
GUT symmetry breaking would be topologically stable 
and accelerated along magnetic field lines. The universe 
is full of long range magnetic fields that would act upon 
the monopoles over their lifetime, likely imparting en- 
ergies ~ 10 15 GeV H. Therefore, magnetic monopoles 
below this energy scale should reach and travel through 
the Earth at relativistic speeds. A relativistic magnetic 
monopole moving through a transparent medium would 
produce copious amounts of Cherenkov light, ^8300 
times a single muon in ice Q. Thus, large Cherenkov 
detectors like IceCube are an ideal experiment to search 
for these particles. 

The current best limits on the flux of magnetic 
monopoles at the 90% confidence level (C.L.) for rel- 
ativistic speeds between /3 = 0.8 and Lorentz boost 
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7 = 10 7 are set by the ANTARES detector @ at the 
~ 10 _17 cm _2 sr _1 s _1 scale. This recent result is the 
first to surpass the results from the AMANDA detec- 
tor Q, IceCube's proof of concept, which set flux limits 
- 3 x KT^cm-^sr-^- 1 . ANTARES also searched for 
magnetic monopoles below the Cherenkov threshold but 
still energetic enough to knock off electrons that pro- 
duce Cherenkov light. This extension sets flux limits at 
the ~ 5 x 10 _17 cm _2 sr _1 s~ 1 scale down to a speed of 
j3 = 0.625. For lower speeds, MACRO provides com- 
prehensive flux limits ~ 10 _16 cm _2 sr _1 s _1 [7] for speeds 
down to f3 ~ 10~ 4 while flux limits at ultra-relativistic 
speeds are set by radio detectors RICE Q and ANITA 
9] at the ~ 10 _19 cm _2 sr _1 s _1 scale. 

These are important as they are flux limits below 
the 'Parker Bound' [10| (~ 10 _15 cm _2 sr _1 s _1 ), an as- 
trophysical flux limit derived by considering the sur- 
vival of the galactic magnetic field in the presence of 
magnetic monopoles. More sophisticated calculations 
that consider velocity [HI relax the bound on rela- 
tivistic magnetic monopoles above a mass of 10 11 GeV 
due to the shortened time spent in the galactic field. 
However, an 'Extended Parker Bound' found by con- 
sidering the survival of a modeled seed field still pro- 
duces flux limits well below experiments, with $ ~ 
10- 16 (Mass M p)/(10 17 GeV) cm^sr^s- 1 

This paper describes the search for relativistic mag- 
netic monopoles in data taken with the IceCube detec- 
tor between May 2007 and April 2008. The analysis is 
optimized for magnetic monopoles with modest Lorentz 
boosts (7 < 10) and charge g — 1. The derived flux limits 
are conservative upper bounds for magnetic monopoles 
with larger 7 or charge, as these cases produce more light 
in the ice. The paper is organized as follows. Section [TT1 
describes the IceCube detector. Section HTT1 describes the 
simulation of background and signal. Section IIVI defines 
the variables and outlines the steps used to discriminate 
signal events from background. Section [V] summarizes 
the uncertainties. Section PVTl presents the final results for 
an isotropic flux of magnetic monopoles at the detector. 
Section IVIII extends this result to an isotropic flux at the 
Earth's surface by considering the energy loss of mag- 
netic monopoles through the Earth. These results are 
presented as a function of both the magnetic monopole 
mass and kinetic energy, allowing for the broad range of 
uncertainty in each of these values. Concluding remarks 
are presented in Section IVIIII 



II. ICECUBE DETECTOR 

IceCube is a telescope at the South Pole which detects 
neutrinos by measuring the Cherenkov light from sec- 
ondary char ged particles produced in neutrino-nucleon 
interactions A total of 5160 Digital Optical Mod- 

ules (DOMs) are arranged in 86 vertical strings frozen in 
the ice between 1500 and 2500 m below the surface over 
a total volume of 1 km 3 . Construction was completed in 



December 2010. The data for this analysis were taken 
during the construction phase, when only 22 of the 86 
strings had been deployed. The 22 strings contain a vol- 
ume of ~ 0.2 km 3 . 

The DOM is the centerpiece of the IceCube detec- 
tor and houses a 10-inch photomultiplier tube (PMT) 
to detect light, onboard electronics for pulse digitiza- 
tion, and LED light sources for calibration. Light sig- 
nals which pass a threshold of 0.25 photo-electron (PE) 
PMT pulse heights are digitized and the DOM is said to 
'launch'. Two types of waveform digitizers are utilized. 
The Analog Transient Waveform Digitizer (ATWD) bins 
the waveform with a 3.3 ns sampling period over a read- 
out window of 420 ns. It supports three channels with 
different gains in order to extend its effective dynamic 
range. The PMT Analog to Digital Converter (ADC) col- 
lects data at a slower sampling rate of 25 ns and records 
for 6.4 fis. 

A time is calculated for the launch by re-syncing the 
threshold crossing to the next leading edge of the inter- 
nal DOM clock, which oscillates at 40 MHz. More pre- 
cise timing is achieved by later reconstructing the leading 
edge of the digitized waveform, though in this analysis 
the coarse time is sufficient. For more information on 
the DOM and its components, see [bH [l5|. 

Each waveform digitizer outputs the signal in terms 
of counts/bin values that directly map to the voltage 
recorded. For the PMT ADC, which is the only chan- 
nel used in this analysis, a single photo-electron corre- 
sponds to ~ 13 counts [3. The PMT ADC saturates at 
1024 counts, which occurs when ~ 50 — 100 PE's arrive 
in a single 25 ns bin. More generally, DOM's that receive 
~ 600 PE's over the full readout window typically satu- 
rate. Figure [T] shows example PMT ADC waveforms for 
both background and signal events at various distances 
to the DOM. The flattened top for the signal indicates 
the point where the digitizer saturates. 

Once a launch is recorded, the DOM checks the four 
nearest neighbors on the string to see if another hit oc- 
curred within a 1/xs time window. By requiring com- 
panion launches to occur, the effect of dark noise hits 
is reduced. If this local coincidence condition is met, 
the digitized waveforms and time are sent to the surface. 
A trigger algorithm is applied to determine if a physics 
event has been detected. For the 22 string detector, this 
algorithm checked if 8 hits were recorded within a sliding 
5 fj,s time window. For data in this analysis, the aver- 
age trigger rate is ~ 550 Hz and is vastly dominated by 
muons generated in cosmic ray air showers in the atmo- 
sphere above the South Pole. 

Against this background a magnetic monopole event 
would stand out due to the much higher light deposi- 
tion. For this analysis a further filter is applied to the 
data online at the South Pole, requiring at least 80 DOM 
launches in an event [llj]. This retains all bright events, 
regardless of direction. The passing rate for this filter is 
~ 1.5 Hz. It consists of muon bundles containing hun- 
dreds of muons generated by high energy cosmic ray pri- 
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FIG. 1: Averaged PMT ADC waveforms for 1000 simulated events of atmospheric muon background (left) as well as the 
f3 = 0.995 (middle) and /3 = 0.8 (right) magnetic monopoles. Distances are how far away the particle is from the DOM 
receiving the light. 



maries. All data that pass this filter are sent north via 
satellite to a data warehouse for use by the entire collab- 
oration. 



III. SIMULATION OF DATASETS 

Simulation of the background and signal are done 
within the ICETRAY framework, a C++ based code written 
for use by the IceCube Collaboration. It includes tools to 
simulate the detector response to light produced by par- 
ticles as well as the triggering and filtering algorithms. 
This allows simulated events to be compared directly to 
the experimental data. 



A. Background Datasets 

Background simulation is composed of muon bundles 
and neutrinos produced in the atmosphere by high energy 
cosmic rays. The generation uses importance sampling 
in energy so that at the final analysis level the statistical 
uncertainty in background prediction is of the order of 
systematic uncertainty or less. 

Atmospheric muon bundles are simulated with 
CORSIKA [17| using two primary types: proton to rep- 
resent light elements and iron to represent heavier 
ions. Primary energies are simulated between 10 4 and 
10 11 GeV. Events are generated with an E~ 2 spectrum 
to oversample the high energy region. The events are 
weighted to fits of extensive air showers introduced by the 
KASCADE Collaboration QJ]. The muon bundles are 
then propagated through the ice using MMC(Muon Monte 
Carlo) pj|. 

ANIS(A11 Neutrino Interaction Simulation) [2(| is used 
to simulate both muon and electron neutrino events. 
The neutrinos are generated with an E~ x spectrum and 
given weights corresponding to a conventional atmo- 
spheric neutrino flux from Honda plj and a prompt flux 
from charmed meson production based on the Enberg, 
Reno, and Sarcevic model 1231. 



B. Signal Datasets 

Code developed specifically for this analysis is used to 
generate and propagate the signal magnetic monopoles. 
Three datasets are created for discrete speeds of /? = 
0.8,0.9 and 0.995(7 — 10)- Monopole tracks are gener- 
ated by randomly distributing vertices on a circular "gen- 
eration plane" with radius 650 m at a distance of 1000 m 
from the detector center. From the vertices, monopoles 
are propagated towards and through the detector with di- 
rections perpendicular to the plane. During generation, 
the orientations of the generation plane relative to the 
detector are randomized, thereby creating an isotropic 
monopole flux through the detector. 

Above j3 ~ 0.1 and below 7 ~ 10 4 , the electromagnetic 
energy loss of magnetic monopoles through matter is well 
described by a combination of ionization and atomic exci- 
tations, collectively referred to as 'collisional' energy loss 
[23l ]. As the choice of simulated events only reach 7=10, 
this is the only energy loss considered in propagation. 
Above 7 ~ 10 4 , energy losses from pair production and 
photo-nuclear interactions surpass the collisional losses. 
These energy losses are considered in Section IVII Al for 
magnetic monopoles traveling through the Earth with 
large boost factors. Bremsstrahlung, which is propor- 
tional to 1/M 2 , is heavily suppressed. 

For each dataset, 100,000 events are generated at a 
mass of M = 10 11 GeV. The effect of choosing one mass 
is mitigated since the Cherenkov light output only de- 
pends on speed which remains essentially constant over 
the 1.2 km path through the detector. 



IV. 



EVENT SELECTION 



The main strategy employed to select relativistic mag- 
netic monopoles is to look for extremely bright events. 
This is measured by counting the number of DOM 
launches which capture a high charge. High charge DOM 
launches are defined as ones that saturate the PMT ADC 
channel. Figure [5] shows the number of these "saturated 



5 




FIG. 2: The number of saturated hits per event for the sim- 
ulated signal and atmospheric muon background (CORSIKA). 
In addition, the full experimental data set is included. 



hits" (NSAT). To visualize the signal event rates, a flux 
of 5 x 10 _17 cm~ 2 sr~ 1 s _1 is used. 

A secondary strategy is to exploit the arrival direc- 
tions of the incoming particle tracks. The dominant back- 
ground of atmospheric muon bundles can only reach the 
detector from above the horizon. This background can be 
suppressed by focusing on events with arrival directions 
below the horizon. 

Event selection consists of three phases. First, a sim- 
ple filter is applied to reduce the data to a manageable 
size. Then, particle tracks are reconstructed and poorly 
reconstructed events are rejected using quality cuts. At 
the final stage, an optimized cut which maximizes the 
Model Rejection Factor (MRF) [24] is found. To reduce 
experimenter bias, the maximized MRF is found using 
simulated background alone. The resulting cut is then 
applied to the experimental data. 

Table Q] displays the final event rates (in events/year) 
for each of the datasets considered at all levels of the 
analysis. 

A. Track Reconstruction 

Since directional information is used mainly to dis- 
tinguish between up and down-going particles, pointing 
accuracy is only of secondary importance. Contrary to 
most IceCube analyses, which use computationally inten- 
sive likelihood methods to reconstruct the particle tracks 
with sub-degree accuracy, a very fast analytic fit proved 
sufficient for this analysis. 

The fundamental piece of datum used by the recon- 
struction is a "hit", which is defined as the location X 
and a time t of a DOM launch. The track direction and 
particle speed are reconstructed by a least-squares fit of 

the observed hit pattern j X- u ti j to a plane wave of light, 

whose analytic solution is given by [25j | 



FIG. 3: Distance from particle track to saturated DOM. All 
histograms are normalized to I. For CORSIKA, a track is de- 
fined by the primary cosmic ray. For data, it is defined by 
the reconstructed track. Events with only two saturated hits 
will reconstruct through the hits and produce the large spike 
at zero. 

X = A avg + Vt (1) 

y _ J2 (jfj ~ A avg )(ti - t avg ) 
(Pi ~ ^avg) 2 

where A avg and t avg are the average position and time 
of all the hits. The hit times ti correspond to the time 
at which the DOM records a launch. Studies of the re- 
construction accuracy demonstrated this to be a better 
definition for ti than the peak time of the PMT pulse, 
likely because the launch time corresponds to the arrival 
time of those Cherenkov photons which are least delayed 
by scattering in the ice. The reconstructed track direc- 
tion is defined by the velocity vector V. 

Because of the simple straight line hypothesis, and be- 
cause the lincfit does not take into account photon prop- 
agation through the ice, the reconstruction accuracy im- 
proves if only hits close to the particle track are included 
in the fit. This is achieved by selecting hits in which a 
large number of photons are detected. A zenith angle 
resolution (~ 2°), defined here as the median difference 
between the true and reconstructed zenith direction for 
simulated events at the penultimate cut level, is achieved 
by only using hits that saturate the PMT ADC. Shown 
in Fig. [3] are the distances from the primary track to a 
saturated hit. Saturated hits are up to ~ 10 m away for 
muons and up to ~ 60 m for the fastest monopoles. The 
relative closeness of the saturated hits mean timing infor- 
mation will be less affected by scattering and absorption, 
improving the accuracy of reconstructing the particle. 

In addition, all hits in which the saturation occurred 
more than 500 ns after the DOM launched are excluded 
from the fit. These are saturated hits where the launch 
time is caused by something other than the saturating 
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particle, e.g. dark noise. This creates errors in the recon- 
struction since the hit information includes a time well 
before the physics event. Roughly 0.05% of the saturate 
hits are removed by this criterion. 

The robustness of the linefit against timing inaccura- 
cies in the hardware was studied by smearing the hit 
times consistent with the frequency of the internal DOM 
clock. This resulted in a negligible change on the recon- 
struction accuracy (< 1%) and final sensitivity (< 1%). 



B. Low level filter 

The low level filter selects events with high Cherenkov 
light yield by requiring at least two of the hits to saturate 
(NSAT>1). This cut reduces background by ~ 99.5% 
and signal by ~ 10 — 15% and ensures that the minimum 
required two hits are available to reconstruct the track 
direction. 




cos 8 



FIG. 4: Final cut on simulated background. This includes 
atmospheric muon bundles and atmospheric neutrinos. Data 
are histogrammed with bin sizes of NSat=l and cos 6=0.022. 



C. Quality Cuts 

Background events which record several saturated hits 
as a result of a bright secondary cascade result in all the 
hits occurring within a small time interval and being lo- 
cated in a relatively small volume. This results in poor 
reconstructions because of the small lever arm to deter- 
mine the overall directionality of the hits. These are re- 
moved by requiring the saturated hits to occur over at 
least 750 ns. This reduces background by ~ 80%. The 
signal is reduced by ~ 30%, but these generally represent 
poor quality events that only saturate one or two strings. 

A second category of mis-reconstructed events are 
caused when multiple muon bundles travel through the 
detector in a single trigger window. These events are 
problematic for this analysis when one saturates a DOM 
in the bottom of the detector before a second saturates a 
DOM near the top, resulting in an up-going reconstruc- 
tion. The large majority are separated enough in time 
so that the speed of the reconstruction (Eq. [5]) connect- 
ing the two events is very low. These are eliminated by 
removing events with a reconstructed speed below 0.2 
m/ns. A second cut requiring the difference in cosQ be- 
tween the reconstruction on all hits versus saturated hits 
to be within 0.6 of 0.0 eliminates these events that hap- 
pen close in time. The combined effect of these cuts is to 
remove ~ 1% of the background and signal. 



D. Final Cut 

The optimized final selection is a piece- wise, linear cut 
on NSAT and cos 8. Figures [4] and [5] show the distribu- 
tion of simulated background and the fastest monopole 
signal in this plane, along with the cut. The background 
is dominated by atmospheric muon bundles, which have 




FIG. 5: Final cut on f3 — 0.995 signal monopoles. Data are 
histogrammed with bin sizes of NSat=l and cos 0=0.022. 



a rate ~5 orders of magnitude larger than the atmo- 
spheric neutrinos. They are essentially all down-going, 
with the more vertical events producing more saturated 
hits. Hence, the final cut is chosen to be on NSAT with 
an angular dependence: constant in the up-going region 
(cos# < 0.0) and linearly increasing in strength in the 
down-going region (cos 8 > 0.0). These two cuts join 
at cos 8 = 0.0. The two numbers that describe this cut 
are the value of the NSAT cut for the up-going region 
('base') and the linear steepness for the down-going re- 
gion ('slope'). The final cut is given by: 



NSAT > 



base if cos 8 <= 

base + slope * cos if cos 8 > 



(3) 



A scan was made through possible values of the base 
from to 25 in increments of one and the slope from 
to 250 in increments of five. For each possible value, the 
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TABLE I: Event rates in events/year for each dataset at all levels of the analysis. Includes simulated signal, background and 
the experimental data. For signal rates, a flux of 5 x 10 _17 cm _2 sr _1 s _1 is assumed. 



Dataset 


Online Filter 


Low Level 


Quality Cuts 


Final 


Experimental Data 


3.15 x 10 Y 


6.55 x 10 b 


1.21 x 10 b 





Corsika Proton 


7.35 x 10 b 


2.65 x 10 b 


2.93 x 10 4 


3.61 x 10~ 4 


Corsika Iron 


5.14 x 10° 


2.20 x 10 s 


6.19 x 10 4 


4.70 x 10~ z 


Atm Conv 


37.9 


26.4 


13.6 


3.45 x W~' 2 


Atm Proinnt 1/ 


4.9 


2.83 


0.334 


4.12 x 10~ z 


Atm Conv v R 


1.4 


0.967 


8.08 x 10~° 


5.52 x 10" b 


Atm Prompt v £ 


2.0 


1.86 


1.43 x 10 -3 


7.39 x 10~ 4 


Bkgrd Total 


1.25 x 10' 


4.85 x 10 b 


9.12 x 10 4 


0.124 


P = 0.995 


100 


89.2 


63.4 


35.6 


/3 = 0.9 


95.3 


84.5 


60.8 


33.4 


P = 0.8 


81.0 


70.1 


46.5 


22.1 




Baseline Cut 



FIG. 6: Scan of Model Rejection Factors for final cut opti- 
mization. The circle corresponds to the minimum value. 

MRF [24[ is found using the event expectation from sim- 
ulation. Figure [6] displays the result of the scan, showing 
the stability of the minimization. The highlighted value 
corresponds to the minimum with a base of 7 and a slope 
of 150. The final cut resulted in a background expecta- 
tion of 0.124 events/year and signal efficiencies ranging 
between ~ 47% to 56% relative to the penultimate cut. 



V. UNCERTAINTIES 

Uncertainties were studied largely with Monte Carlo 
simulations. Table |TT] contains the results. The large 
relative background uncertainty is acceptable given the 
small absolute event rate. Uncertainties consisted of 
three types: (1) Theoretical uncertainties in the simu- 
lated models, (2) Uncertainties in the detector response, 
and (3) Statistical uncertainties. 

Theoretical uncertainties include the shape and nor- 
malization of the background energy spectrum for both 



the atmospheric muons and neutrinos. In addition, the 
cross section uncertainty modeled in both MMC and ANIS 
is studied. Detector uncertainties include uncertainties 
in the scattering and absorption parameters of the ice 
and the efficiency of the DOM. 

For atmospheric muon background, the dominant un- 
certainty is from the cosmic ray energy spectrum. For 
both elements, the parameters of the assumed broken 
power-law (break energy, power-law indices below and 
above the break, and absolute normalization) were var- 
ied within the uncertainties in the 2-component model 
18]. For iron, the extreme case of no break is taken as 
the upper end of the uncertainty since the expected break 
occurs beyond the fit region of the model. Since the fi- 
nal CORSIKA sample is overwhelmingly high energy iron 
primaries above 10 10 GeV, it is very sensitive to changes 
in the spectral weighting values. 

The conservative nature of this assumption is to al- 
low for uncertainties at the high energy range that are 
not easily tested by simulation. Despite this extreme, 
the absolute uncertainty is still less than 0.5 events/year. 
Signal is more robust due to the brighter light yield. This 
allows a larger sample to pass the final cuts relative to 
background causing it to be less sensitive to variations in 
the detector response. 



VI. RESULTS 

The optimized cut is then applied to the full experi- 
mental data sample. No events survived on an expected 
background of 0.124 events, resulting in a Feldman and 
Cousins upper limit of 2.44 at the 90% C.L. [26]. The 
final distribution is shown in Fig. [7] 

The final flux limit is calculated incorporating the sys- 
tematic and statistical uncertainties using the profile log- 
likelihood method implemented in the P0LE++ program 

Table Mil displays the resulting sensitivities and final 
limits on the flux of magnetic monopoles at the detector 
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TABLE II: Relative uncertainties for predicted event rates of background and signal. Total uncertainties found by adding 
absolute rate deviations in quadrature. 



Uncertainty 




Backf 


,'round 




Signal 




CORSIKA 




Ve 


Total 


P = 0.8 


P = 0.9 


P = 0.995 


Normalization 


26% 


1 1 (V 

11% 


< I /O 


12% 




Spectrum 


990% 


22% 


39% 


380% 




MMC cross-section 


10% 


10% 




7.4% 




v cross-section 




6.4% 


6.4% 


4.0% 








DOM Efficiency 


27% 


38% 


38% 


25% 


5.8% 


4.8% 


1.4% 


Ice Properties 


78% 


40% 


71% 


40% 


7.1% 


4.2% 


0.2% 


Statistical 


22% 


14% 


19% 


12% 


0.9% 


0.7% 


0.7% 


TOTAL 


990% 


64% 


110% 


382% 


9.2% 


6.5% 


1.7% 




1(T 



10" 



1(T 



10" 





pown 




AMANDA 




\up — 


_ _ JVIACRQ 






^ ANTARES 

-_____lceCube 




RICE/ANITA^ 



0.6 



0.7 



0.8 



0.9 



1.0 



FIG. 7: Final cut on the full experimental data sample. Data 
is histogrammed with bin sizes of NSat=l and cos 0=0.022. 



FIG. 8: Flux limits (90% C.L.) at the detector as a func- 
tion of p. Includes results from AMANDA ((J, Baikal HH, 
ANTARES [I, MACRO 0, RICE @], and ANITA [|] 



TABLE III: Final sensitivities and limits (90% C.L.) on the 
flux of magnetic monopoles at detector in cm _2 sr _1 s _1 







P = 0.8 


P = 0.9 


P = 0.995 


Sensitivity 


6.10 x 10" i5 


3.94 x 10" is 


3.73 x 10" ia 


Final Limit 


5.57 x 10~ is 


3.56 x 10 _i8 


3.38 x 10 _1S 



at the 90% C.L. Figure [5] shows this result compared with 
previous searches from neutrino telescopes. 



VII. DISCUSSION 

In order to describe the results as pertaining to an 
isotropic flux at the surface of the Earth, the efficiency of 
the analysis as a function of zenith is combined with the 
acceptance of rclativistic magnetic monopoles through 
the Earth. The previous AMANDA analysis did a similar 
procedure @. 



A. Angular Acceptance Through the Earth 

For an isotropic, mono-energetic flux $T s ' M of mag- 
netic monopoles with mass M and kinetic energy E^in = 
M(-f s — 1) at the Earth's surface, the resulting 7 of the 
monopole flux at the detector is calculated for cos 9 val- 
ues in increments of 0.1. The energy loss is modeled 
using Ahlen's stopping power formula for collisional loss 
[HI and code from ANITA Q for pair production and 
photo- nuclear losses. 

Figure [5] shows the angular acceptance of relativistic 
magnetic monopoles traveling through the Earth. Each 
line indicates the boundary between the mass and kinetic 
energy values which allow the monopole to reach IceCube 
at a particular speed threshold for a given zenith. For in- 
stance, if the mass and kinetic energy are in the region 
above the cos6> = —1.0 line, these describe a magnetic 
monopole that can remain relativistic traveling the di- 
ameter of the Earth, while those above the cos 9 = 1.0 
can remain relativistic traveling through the atmosphere 
and the ~ 2 km of South Polar ice to reach the detector. 

The shape of the lines can be understood by consider- 
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Monopole Mass (GeV) 



FIG. 9: Parameter space of magnetic monopoles at the 
Earth's surface required to reach the IceCube detector at rel- 
ativistic speeds. Dotted is for a speed threshold of ft = 0.995, 
dashed for /3 = 0.9, and solid for /3 = 0.8. 



ing the full acceptance (cos 9 = —1.0) case: 

* The collisional energy loss straight up through the 
Earth is ~ 10 11 GeV. This loss is not enough to slow 
rclativistic magnetic monopoles with masses above 
~ 10 12 GeV to sub-relativistic speeds. Therefore, 
the acceptance is determined solely by the starting 
energy. 

* Magnetic monopoles with masses between ~ 10 7 — 
10 12 GeV can still reach the detector so long as 
there is enough kinetic energy to overcome the col- 
lisional loss. Hence, the line flattens out around 
~ 10 11 GeV. 

* For magnetic monopoles with masses below ~ 
10 7 GeV, the necessary starting energy begins to 
increase to overcome the increasing effect of pair 
production and photo-nuclear energy losses, which 
begin to dominate for 7 ~ 10 4 . 



B. Angular Acceptance of Analysis 

The analysis is much more sensitive to an up-going 
signal, due to the large atmospheric muon bundle back- 
ground. This is described quantitatively by calculating 
the effective area as a function of zenith. The effective 
area corresponds to the cross sectional area of an ideal 
detector with 100% efficiency. Using the same cos 9 bins 
as above, the effective area is given by: 



A^ ff (cos( 



Njjpos 

gen tvt7 

1 Vp-f 



{cos9) 



(4) 



where A^ cn is the area of the generation plane for a given 
7 and Nj et /N^ en is the fraction of magnetic monopoles 
generated that survive the final analysis cut. 



FIG. 10: Effective Area for each cos 6 bin. 



Figure [TU] shows the result of the three generated 
speeds. From Section IIIIB| AJ en = 1.33 km 2 and 
7V^ en (cos 9) = N^ en /20 = 5000, since the generated flux is 
isotropic. This is conservatively generalized to any speed 
at the detector by treating the effective area as a step 



function, 



e.g. 



A eff 



effective area is set to zero. 



17=10 

off 



etc. For f3 < 0.8, the 



C. Limits on Isotropic Fluxes at the Earth's 
Surface 



The final limit on a flux with given mass and kinetic 
energy at the Earth's surface ($2o%C l ) * s calculated by 
scaling a reference flux with the ratio of the Feldman- 
Cousins upper limit (^90%) [IB] to the expected number 
of signal events seen in the detector using the reference 
flux. The expected signal event number is found by going 
through each cos 9 bin and determining (1) what speed 
the monopole will have at the detector (7^) based on 
Section IVII Al and (2) calculating the effective area for 
that speed and cos 9 bin based on Section IVII Bl The 
final flux limit becomes: 



4> 



7s, M _ 
90%C.L. ~~ 



v si 



/x 90% (V obs = 0.0,V, 



bkg 



0-124) 7.,a/ 

^Rcf 



(5) 



20 

3 $£ e f 2tt £ (A cos 9)iA& (cos Bi) (6) 
i=l 



Xii vc = 2.06 x 10 7 s is the total livetime of the analyzed 
data set, N^ g = 0.124 is the final background expecta- 
tion from Table HI A cos (9 = 0.1 is the width of the cos# 
bins, and the 2-7T arises from the azimuthal symmetry of 
the Earth. For most tested values of j S) M, the final 
speed is the same for all bins and the flux limit calcula- 
tion returns the same answer as Table IIII1 

To place this result in context, Fig. [TT] displays the 
current best experimental flux limits over a wide range 
of mass and kinetic energy values of magnetic monopoles. 
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log 10 (M M oNOPOLB) GeV 

FIG. 11: Final flux limits (90% C.L.) as function of monopole 
mass and kinetic energy at Earth's surface. For relativistic 
mass and energies, only the most restrictive limit is displayed. 
Includes the Parker Bound [10J and results from MACRO [7(, 
RICE [3], and ANITA [f|. For numerical values of the final 
result for this analysis, see [2!|. 



Below 7 = 1.67 the analysis does not apply as the 
monopoles fall below the Cherenkov threshold, while 
above 7 = 10 7 , the radio neutrino detectors offer bet- 
ter sensitivity. 

For the range of mass/kinetic energy pairs resulting 
in 1.67 < 7 < 10 7 , this analysis provides in general the 
best flux limits to date. The exception occurs for the 
smallest masses and kinetic energies, where attenuation 
in the Earth affects the signal acceptance. To help guide 
the eye, lines showing the angular acceptance solid an- 
gle il for the (3 = 0.9 magnetic monopoles are included. 
The solid angle is found by multiplying 2-k by the range 
of cos 9 for which the mass and energy combination can 
reach the detector. Hence the shape matches Fig. [9] 
As the solid angle approaches 2n, acceptance below the 
horizon is lost and the limit becomes much weaker. 

For the cases where 7 > 10 4 , the flux limit from this 
analysis is conservative, as the monopole would have a 
large light contribution from secondary cascades which 
are not yet included in the simulation. These will make 
the event brighter in the detector and increase the selec- 
tion efficiency. 



VIII. CONCLUSION 

This analysis is the first search for magnetic monopoles 
using the next generation of neutrino telescopes. A final 
flux limit of $ 90 %c.L. = 3.38 x lO-^cm^si'-V 1 for /3 > 
0.995 is found. For speeds down to j3 = 0.8, the flux limit 
is slightly higher. This applies to an isotropic flux at the 
Earth's surface for relativistic magnetic monopoles with 
mass above ~ 10 6 GeV and energy above ~ 10 10 GeV 
(Fig. ITT]) . Even with a single year of data operating at 



~ 20 % of the final instrumented volume, experimental 
flux limits are achieved that are a factor of 4 below the 
current best constraints up to 7 ~ 10 7 and provide a 
good compliment to the more sensitive radio searches for 
ultra-relativistic monopoles. 

This analysis does not follow IceCube's usual proce- 
dure of a blind analysis. An original analysis performed 
on the data was done in a blind fashion, with cuts be- 
ing determined by simulation datasets along with a 10% 
'burn' sample of experimental data. It aimed to enhance 
the sensitivity to slower monopoles by binning the data 
based on speed reconstruction. Unblinding revealed de- 
ficiencies in the background simulation to reproduce the 
tails of the speed distribution where the slower signal 
should be and allowed obvious background events into 
the final sample. After determining no monopole events 
were recorded, the analysis reported here is performed 
with cuts optimized on improved simulation and not the 
experimental data. The only changes involve a slight 
tightening of quality cuts motivated by the new simula- 
tion and abandoning the binning based on speed recon- 
struction. For a full description of the original analysis, 
final event rejection, and motivation for changes, see [29j. 

Preliminary work on the 2008 data run shows that the 
increased detector size and improvements to the analysis 
method provide a further factor of 4 reduction in the 
sensitivity [30( | . With more data and refined techniques, 
IceCube and other neutrino telescopes will continue to 
prove valuable in searches for magnetic monopoles in the 
relativistic regime. 
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